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The ground-state phase diagram is numerically studied for an electronic model consisting of the 
spin exchange term (J) and the correlated hopping term (J3: the three-site term). This model has 
no single-particle hopping and the ratio of the two terms is controlled by a parameter a = 4t3 /J. 
The case of a = 1 corresponds to complete suppression of single-particle hopping in the strong- 
coupling limit of the Hubbard model. In one dimension, phase separation takes place below a 
critical value a c = 0.36-0.63 which depends on the electron density. Spin gap opens in the whole 
region except the phase-separated one. For and low hole densities, charge-density-wave 

correlations are the most dominant, whereas singlet-pairing correlations are the most dominant 
in the remaining region. The possibility of superconductivity in the two-dimensional case is also 
discussed, based on equal-time pairing correlations. 
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The discovery of high-T c cupratesB has driven us 
the question whether the new microscopic theory for 
pairing mechanism is required. In the cuprates strong 
Coulomb repulsion may be essential, and various pair- 
ing mechanisms have been proposed under strong cor- 
relation. One of the common features in experiment 
is the extremely flat band around the momenta (0, ir) 
and (n, 0) in units of the reciprocal lattice parameter, 
which has been observed-in angle-resolved photoemission 
spectroscopy (ARPES).B This feature is attributed to a 
strongly correlated effect. The band structure similar to 
the ARPES data has been numerically obtained for the 
Hubbard model with strong on-site Coulomb repulsion© 
and the t-J modeler in two dimensions. 

If the flat band plays a crucial role in increasing the 
critical temperature, proper control of the single-particle 
process should lead t^more enhancement of supercon- 
ducting correlations .o&Lr Numerical results for one- 
dimensional (ID) electron systems show that flattening 
of the band near the Fermi points actually extends the 
parameteij-region where pairing correlations are the most 
dominant \y If single-particle hopping is completely ab- 
sent as inessentials, other hopping mechanism is neces- 
sary for superconductivity. One of the candidates is the 
correlated hopping, which is often called the three-site 
term and is naturally derived in,thc strong-coupling ex- 
pansion of the Hubbard modelBEP However, in one di- 
mension, the correlated hopping is known to be suppress 
pairing correlations in the low-hole-density region J13'tiP 
This is because the-Jiorrelated hopping leads to a repul- 
sion between holes.td) In contrast, spin exchange inter- 
action favors the attractive configuration of electrons. 
Therefore the correlated hopping and the spin exchange 
should compete in an electron-density region, and the 
ground-state phase diagram is nontrivial. 
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In this paper, we shall clarify the ground-state prop- 
erties due to the competition between spin exchange and 
correlated hopping in one dimension. We also discuss re- 
alization of superconductivity under complete suppres- 
sion of single-particle hopping in the 2D case. To this 
end, exact diagonalization is employed for finite systems. 

In the Hilbert space without double occupancy of elec- 
trons, the strong-coupling expansion of the Hubbard 
model yields the following Hamiltonian up to tiy$ order 
t 2 /U (U being the on-site Coulomb repulsion): 



JL = Jit + JLj + JL 
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where (ij) and (ij£) are nearest neighbors, J = 4t 2 /U 
and ts = a J/4 with a — 1. The constrained fcrmion 
operator is given by Cj CT = Ci CT (l — nj,_o-). Here we 
consider the value of a as a continuous parameter. In one 
dimension, the ground-state phase diagram in the J/t-n 
plane (n beingj-the electron density) has been, investi- 
gated for a = ut3 and several values of a > 0£3 We re- 
fer to a model with Tit deleted, namely Tijtz = T~tj+Ht3, 
as the J-ts model. We will complete a phase diagram 
in the a-n plane for the J-ts model. Originally, both 
Hj and Ht3 decay more rapidly than Jit when the hop- 
ping integral t approaches zero. In considering possible 
mechanisms for superconductivity, however, it would be 
meaningful to deal with each term of eq. (1) individually 
away from the Hubbard model. 
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It should be noted that in the ID J-ts model the set 
of all eigenvalues under the antiperiodic boundary con- 
dition is equivalent to that under the periodic one. 

We first discuss the phase diagram in the ID case. In 
order to estimate the boundary of phase separation, we 
calculate the compressibility given by 

N 4 
K ~ N* ' E (N; N e + 2) + E (N; N c - 2) - 2E (N- N e ) ' 

(5) 

where Eo(N;N e ) denotes the ground-state energy with 
iV e electrons in N sites. In Fig. 1 we show the values 
of a c at which 1/ n crosses zero for TV = 10-16. The size 
dependence is small and this implies a reliable estimate of 
the thermodynamic results. The phase separation occurs 
for a < a c . 

In the low-electron-density limit, one can estimate the 
phase boundary by solving two-electron problem.tir In 
the J-t 3 model, the phase separation disappears at n ~ 
for a > a c = 2 In 2 - 1 ~ 0.386© This estimate would 
be exact if there are no bound states with more than 
two electrons. In the high-electron-density limit, on the 
other hand, one may estimate the boundary from bal- 
ance between the gain in kinetic energy and the loss of 
exchange energy due to insertion of two holes in a Heisen- 
berg chain.O Following ref. ,113 we calculate the ground- 
state energy with N — 2 electrons in N sites relative to 
the Heisenberg energy in N — 2 sites. This energy be- 
comes negative for a > 0.497, 0.488, 0.483 and 0.480 for 
N = 10, 12, 14 and 16, respectively. We obtain a c ~ 0.47 
from extrapolation to the thermodynamic limit. This 
value is smaller than a c ~ 0.63 obtained from the com- 
pressibility at the highest electron density (n = 0.875) in 
our calculation. Anyway the critical values at n ~ and 
n ~ 1 are both smaller than a — 1, which implies that 
the a = 1 line is not interrupted by phase separation. 
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Fig. 1. Phase diagram for the ground state of the ID J-t^ model. 
PS: phase-separated state; CDW: charge-density-wave state; SP: 
singlet-pairing state. Closed symbols and open ones show the 
phase-separated boundary and the contour of K p = 1, respec- 
tively, for N = 10 (upward triangles), 12 (squares), 14 (down- 
ward triangles) and 16 (diamonds with a solid line). The con- 
tours of other values of K p for N = 16 are shown by solid 
lines. The low-density limit of the phase-separated boundary! 1 y 
is shown by a broken line. 



We next calculate the spin gap given by Esq = 



limjv^oo-EGs(iV;iVe = Nn), where E SG {N;N C ) = 
E (N;N e ;S? ot = 1) - E (N; N e ; S t ^ = 0). Here S z tot 
denotes the z-component of total spin for the system. 
We estimate the magnitude of spin gap with use of the 
data of N = 9,12,15 for n = 2/3, N = 8,12,16 for 
n = 0.5, and N = 6, 12, 18 for n = 1/3. The following 
fitting function is used for extrapolation to N — > oox3 

T(N) =T{oo) + A/N m + B/N 2m , (6) 

with r = Eqq and m = 1. Figure 2 shows the extrapo- 
lated values of the spin gap as a function of a. At any 
electron density, the presence of spin gap is clearly shown 
for the whole region of a > a c . In the inset of Fig. 2, 
both horizontal and vertical axes are scaled by a. This 
figure makes the continuity to a case of only TL t z (i.e., 
1/a = 0)11^ clear. The spin gap scaled by a has a peak 
near a — 1 at any filling. 
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Fig. 2. Spin gap versus a for the ID J-ts model. Each value 
is extrapolated from the finite-size data. The values for n = 
2/3,1/2 and 1/3 are shown by upward triangles, squares and 
downward triangles, respectively. The inset shows the spin gap 
scaled by a as a function of 1/a. The values— «f the analytic 
result for only the correlated hopping (1/a = OH 1 1 P are indicated 
by closed symbols. 

An accurate way to estimate a spiri-gap region is the 
singlet-triplet level crossing methodO'EEP Nakamura et 
al. applied the method to the ID t-J model, assuming 
that the low-energy behavior of the model is described 
by the effective Hamiltonian consisting the U(ll Gaus- 
sian model and the SU(2) sine-Gordon modelJiS One 
can check the validity by seeing that x r = (x ss + 3x s t)/4 
is close to 1/2, where x ss (x st ) denotes the scaling di- 
mension of the singlet (triplet) excitations. These ex- 
citations are calculated under the boundary condition 
twisted from the one for the ground state.EJf 

When the level crossing method is applied to the J-t% 
model, we should notice two points. One is the presence 
of correlated hopping, and the other is the absence of 
single-particle hopping. Here we assume that in the low- 
energy region the model given by eq. (1) is described by 
the above effective Hamiltonian even for t% > 0. The va- 
lidity is checked from the scaling dimension to be shown 
later. Meanwhile, the presence of single-particle hop- 
ping t seems to be indispensable because it provides the 
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finite gradient for linearization of the band dispersion. 
To avoid the ambiguity, we begin with the finite value of 
t, and then consider the limit of t — ► 0. 

Figure 3(a) shows the spin-gap boundary obtained by 
the level crossing method, in the a-t/J plane at n = 0.5 
and 0.875 in 16 sites. There is no spin gap in the region 
above the solid line [case (A)], while the remaining region 
should be the spin-gap region or the phase-separated one 
[case (B)]. Obviously the limit of t — > belongs to the 
case (B) at both electron densities. The validity of the 
method is ensured in Fig. 3(b), which shows that the 
scaling dimension x r in the gapless region (t/J^0.8) is 
very close to 1/2 (0.502-0.504) for N = 16,n = 0.5 and 
a = 1. One may ask whether the limit of t — ► and 
the case of t = are connected smoothly. In Fig. 3(c) 
we show the magnitude of spin gap as a function of t/ J. 
The figure informs us no singularity of the spin gap at 
t/J = 0. This situation would hold also for other values 
of a and fillings unless phase separation occurs. Thus 
we conclude that the J-t% model has a spinjjap in the 
whole region except the phase separated oneJiZP 
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Fig. 3. (a) Spin-gap boundary in the a-t/J plane at n = 0.5 and 
0.875. (b) Scaling dimension of the lowest singlet excitation (x SB ) 
and that of the triplet excitation (a; s t) for TV = 16, n = 0.5 and 
a = 1. The quantities x ss , x s t and x r = (x BB +3x a t)/4 are shown 
by circles, squares and diamonds, respectively, (c) Magnitude of 
spin gap as a function of t/J for n = 0.5 and a = 1. Each value 
is extrapolated from the finite-size data (TV = 8, 12, 16). 



Although the J-ta model has no single-particle hop- 
ping, the charge velocity is not zero but finite due to 
the correlated hoppingO Thus we expect that charge 
degrees of freedom for the system behave as a Tomonaga- 
Luttinger (TL) liquid. To support this, we calculate the 
central charge c which characterizes the universality class 
of a model. One can obtain c by fitting the ground-state 
energy as eq. (^) with T = Eq/N and m — 2. Here 
the coefficient A is given by —nv c c/6, and the charge 
velocity v c is extrapolated by eq. (^) with T — v c and 
m = 2. Note that A does not include the spin velocity 
due to the presence of spin gap. With use of the data of 
N = 8, 12 and 16 with n = 0.5, we find that c is almost 
unity (c ~ 0.999) in the region a — 0.6-4.0. This in- 



dicates that charge degrees of freedom are described by 
a TL liquid. For aSaO.4, c deviates largely from unity, 
reflecting the phase separation which is no longer a TL 
liquid. 

In the region c = 1, one can use the following universal 
formula about the correlation exponent K p : 



K p = n 



Dn 2 n 



where the Drude weight D is defined as 
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Here Eo(cf>) denotes the ground-state energy of the sys- 
tem under the twisted boundary condition (ci + N_ a — 
e'^Cia). When a spin gap is present, singlet-pairing 
[charge-density-wave (CDW)] correlation functions de- 
pend on the distance r between sites as r~ 1 / Kp [r~ K >'] 
for r 3> 1. If K p > 1, singlet-pairing correlations are 
the most dominant; otherwise CDW correlations are the 
most dominant. In Fig. 1 we show the contour lines of 
K p with 1.0-2.0. 

Our result for the values of K p is supported by direct 
calculation of the equal-time pairing correlation function 
C(r) = {l/N)Y, t {P}Pi+r)- Here the singlet pairing op- 
erator P t is given by Pi = (0^+14 - c a cj + i iT )/V2, 
and denotes the expectation value in the ground 

state.ES' Figures 4(a) and 4(b) show the distance de- 
pendence of C(r) for various values of a with n = 0.75 
fixed and that for various fillings with a = 00 fixed, re- 
spectively. In any case, pairing correlations reach up to 
long distance for K p > 1. The sum S = ^ r C(r) pro- 
vides the uniform (zero-momentum) component of the 
pairing structure factor. Figure 4(c) shows the result of 
5 as a function of a for various fillings. The values of 
a at which S has a peak are found to correspond to the 
phase-separated boundary (K p = 00) in Fig. 1. 

The phase diagram consists of three regions: the 
phase-separated region, the CDW one and the singlet- 
pairing one. Note that for a — 00 the CDW correla- 
tions ace the most dominant in the range n > 2 — \J~2~ ~ 
0.586.liiP We emphasize that the pairing correlations are 
the most dominant irrespective of n for a = 1 corre- 
sponding to complete suppression of single-particle hop- 
ping in the strong coupling limit of the Hubbard model. 

How is the situation changed on a square lattice? We 
employ a 4x4 cluster with the periodic boundary condi- 
tion in both x and y directions. To search indications of 
superconductivity, we calculate equal-time pairing cor- 
relations given by C s(d) (f) = (1/N) ^ ? (A t ± (?)A ± (? + 
f))BH The pairing operator A±(i) is defined as 
A±u) = &v(o, ? 1 +Qr ? licy,? ±0 ? ,), where + 

= V I %\ \ l+t> X ,l I— t) X ,l l+Oy,l l—Oy,l" 

(— ) corresponds to extended-s-wave (d^^^-wave) sym- 
metry, and the unit vector in x (y) direction is repre- 
sented by 4 (S y )B> In Fig. 5(a) we show the distance 
dependence of C s (r) and Cd(r) with r = \f] for n = 0.5 
and a = 1. The extended-s-wave pairing correlations axe 
stronger than the d;,,2_ a 2-wave ones at all distances. EP 
Figure 5(b) shows the extended-s-wave correlations for 
n = 0.5 and various values of a. The rapid decay at 
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Fig. 4. (a) Equal-time pairing correlation function as a function 
of distance for the ID J-tz model with N = 16, n = 0.75 and 
three different values of a. (b) Same as (a) but with a = co and 
three different values of filling. In (a) and (b), cases indicated by 
closed symbols are within the singlet-pairing region (1 < K„ < 
co) in Fig. 1. (c) Uniform component of the pairing structure 
factor as a function of a at n = 0.25 (squares), 0.5 (circles) and 
0.75 (upward triangles) in a 16-site cluster. 



Fig. 5. (a) Equal-time pairing correlation function as a function 
of distance for the 2D J-ts model with N = 4 X 4, n = 0.5 and 
a = 1. Open (closed) symbols denote extended-s-wave (d x 2_ y 2- 
wave) correlations throughout Fig. 5. (b) Extended-s-wave cor- 
relations with three different values of a. Other parameters are 
the same as in (a), (c) Uniform component of the pairing struc- 
ture factor as a function of a at n = 0.25 (squares), 0.5 (circles) 
and 0.75 (upward triangles). 



long distances for a = 0.2 reflects phase separation. The 
behavior for a = 1 and a = oo (i.e., only the corre- 
lated hopping) is almost the same. This means that the 
pairing correlations are hardly affected by the spin ex- 
change unless phase separation occurs. Finally we show 
the uniform component of the pairing structure factor 
Ss(d) = ErC s(d) (r) in Fig. 5(c). For n = 0.25 and 0.5, 
we observe the abrupt increase in S s as a function of 
a. This implies a transition from phase separation to 
s-wave superconductivity. At fixed a(^0.4), as the elec- 
tron density increases the extended-s-wave correlations 
become weak. Instead the d^a^a-wave correlations be- 
come strong, but these do not excel the s-wave ones at 
least for n < 0.75. 

In summary, we have found that the phase diagram 
of the ID J-tz model is composed of the three different 
regions. Most of the parameter space is occupied by the 
singlet-pairing region. Also in two dimensions, extended- 
s-wave superconducting states are likely to cover a wide 
range, although we have not confirmed whether the long- 
range order exists. A characteristic feature of the su- 
perconductivity discussed in this paper is to be caused 
by only the two-particle process. This mechanism is 
essentially different, from the Bardeen-Cooper-Schrieffer 
(BCS) mechanismEJ which is based on the single-particle 
process and an attractive potential between two elec- 
trons. 
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